Scattering mechanisms and Boltzmann transport in graphene 
Shaffique Adam & ^ , E. H. Hwang a and S. Das Sarma a 

'■Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, MD 20742-4111, USA 



Abstract 

Different scattering mechanisms in graphene are explored and conductivity is calculated within the Boltzmann transport 
theory. We provide results for short-range scattering using the Random Phase Approximation for electron screening, as well 
as analytical expressions for the dependence of conductivity on the dielectric constant of the substrate. We further examine 
the effect of ripples on the transport using a surface roughness model developed for semiconductor heterostructures. We 
find that close to the Dirac point, a ~ n@ , where (3 = 1, 0, —2 for Coulomb, short-range and surface roughness respectively; 
implying that Coulomb scattering dominates over both short-range and surface roughness scattering at low density. 
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1. Introduction 

Over the last two years, there has been a proliferation 
in theoretical and experimental interest in graphene. 
Graphene is a two dimensional sheet of Carbon atoms 
arranged in a honey-comb lattice, or alternatively, a 
mono-atomic sheet of carbon that has been cleaved 
from bulk graphite onto a semiconductor substrate. 
Graphene has generated much excitement as theorists 
and experimentalists explore both the familiar and 
peculiar features of this newly discovered condensed 
matter system. Perhaps the most curious property of 
graphene was the "minimum conductivity" puzzle, 
where a naive examination of experimental results 
would suggest that graphene had a finite conductivity 
even as the density of carriers vanishes. This "carrier- 
free" universal minimum conductivity has been the 
subject of much speculation in the recent literature. 



However, there is now mounting evidence |l|; 0) sug- 
gesting a more sensible alternative - charged impuri- 
ties that are invariably present in the substrate create 
an inhomogeneous potential landscape (0). At large 



external gate voltage V g = n/a 3> 



/a, where 
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riimp is the density of charged substrate impurities 
and a is related to the capacitance of the substrate, 
the potential fluctuations can be ignored, but at low 
density, they locally dope the graphene sheet breaking 
the system into puddles of electrons and holes. This 
spatial inhomogeneity has been directly observed in 
recent experiments We observe that the charged 
impurities therefore have a dual role: first, the induced 
graphene carrier density needs to be determined by 
the screened charged impurity potential, and second, 
the conductivity is determined by charged impurity 
scattering. 

Unlike gapped 2D systems - where much of the same 
physics is at play -, in graphene the electrons and holes 
are both conducting. So whereas in semiconductor sys- 
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tems one might have a percolation transition between 
a conducting and insulating state (0) , graphene trans- 
port, through this percolating network of electron and 
hole puddles, is a transition between two conducting 
states. One also finds that the system self-averages at 
length scales of the correlation length Lj, where <C 
L and L ~ 1 /im is the typical system size. This then 
allows us to calculate the transport properties using 
the rms carrier density n* induced by the potential 
fluctuations (for details, see Ref.jf])). This mean-field 
theory is valid so long as n* is sufficiently large so that 
the percolation transport is not critical - a condition 
we believe holds true for experiments currently being 
done on bulk graphene samples. Ref.(0) recently calcu- 
lated the modification of the critical exponents due to 
the tunneling between the electron and hole puddles. 

In the present work, we examine two other types of 
scattering that may be present in graphene, namely, 
short-range scattering (arising from defects or disloca- 
tions in the Carbon lattice) and scattering from ripples. 
Both of these possible additional sources of scattering 
have been observed in recent surface probe measure- 
ments of graphene (Q). These experiments, however, 
also showed that graphene sheets are relatively defect 
free and the ripples, although present, are very small, 
where the average ripple height A = 1.9 A, while the 
in-plane correlation length £ = 320 A. These numbers 
cast serious doubt over whether the electron scattering 
off graphene ripples could have any observable effect 
on the transport properties of graphene. We observe 
that diffraction experiments performed on suspended 
graphene sheets show ripples of about the same height 
but with an order of magnitude smaller in-plane corre- 
lation length (9) further suggesting that the ripples of 
graphene on a substrate is largely determined by the 
roughness of that substrate. 

For our purposes we reserve the term "long- 
range" scattering for Coulomb impurities and use the 
two phrases interchangeably. We consider all delta- 
correlated disorder (i.e. Gaussian white noise poten- 
tials) to be "short-range" scattering. Other definitions 
have been employed in the literature where white- 
noise disorder is sometimes described as long-range if 
it preserves the graphene valley degeneracy and short- 
range only if the disorder breaks this symmetry. In all 
of our discussion below we shall assume that both spin 
and valley degeneracy is always preserved. 

Short-range scattering has been considered previ- 



ously by many groups including Refs. 

(0; El El) 

us- 
ing various approximation schemes. Here, we present 
results for the screened short-range scattering using 
the Random Phase Approximation (RPA) for electron 
screening iflih which allows us to extract the depen- 
dence of the conductivity on the interaction parameter 
r s , which can be tuned by using substrates of different 
dielectric constants. 

Coulomb scattering was considered by Ref. |12j) in 
the completely screened approximation valid only in 
the limit r s — ► oo and numerically in Ref. ifiol : 1 141: 13). 
To compare the different scattering mechanisms with 
each other and to provide the framework for our calcu- 
lation, we reproduce here the analytic Coulomb scat- 
tering conductivity result that was previously reported 
in Ref. (@). 

The scattering of electrons off ripples has been of re- 
cent theoretical interest. Despite the magnitude of the 
ripples being quite small, it is still believed to be re- 
sponsible for the suppression of weak localization jlih . 
as well as causing changes to the local chemical po- 
tential iflrih and modification of the transport proper- 
ties H3). To provide an estimate for the contributions 
of static ripples to graphene transport properties, we 
use a surface roughness model that was developed to 
model interface potential inhomogeneities in semicon- 
ductor heterojunction systems. We do not attempt here 
to provide a rigorous microscopic justification for ap- 
plying this surface roughness model fl8|) to graphene, 
although given its success in modeling surface inhomo- 
geneities in Si-MOSFETs, it provides a good starting 
point to model the matrix element for electrons scat- 
tering off the potential inhomogeneity caused by rip- 
ples. Whereas the formal validity of this model has yet 
to be established, we are persuaded that it captures at 
least qualitatively, many of the features of the effect of 
ripples on graphene transport. 



2. Formalism 

We calculate the Drude-Boltzmann conductivity 
(see Ref. (@) for a detailed discussion on the range of 
validity of this formalism) using a = (2e 2 /h) E f t, 
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where E-p is the Fermi energy, ft _1 7 is the Fermi veloc- 
ity and the momentum transfer q = |k — k'| = 
The different types of impurities change the scattering 
potential v(r]) and screening is accounted for using the 
RPA dielectric function e(rj) reported in Ref. (flih. For 
Coulomb scatterers, we use 
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where Fi(x) is given analytically by 



Fi(x) ir 3x 2 tv 2 arccos[l/s] . , 

— = — -2)- ? _ r . (3) 

It is straight forward to generalize this formalism 
to calculate the conductivity of screened short-range 
scatterers. We have 
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This result demonstrates that screening does not 
change the carrier density independence of the conduc- 
tivity arising from short-range scattering and further 
predicts how the conductivity depends on the dielec- 
tric constant of the substrate. 

In similar fashion, we can solve the surface roughness 
model lflj[) . First to determine its importance near the 
Dirac point, we calculate the limit n — > and find 



a H (n -> 0) = 



1 
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Therefore, within this surface roughness model ( and 
in a manner qualitatively similar to short-range scat- 
tering), this effect becomes unimportant at the lowest 
carrier densities and at low temperature. Since exper- 
imentally, the parameter \/im*^ > 3, the limit 3> 
1 is more realistic. An analytic result can be obtained 
in the unscreened limit (where we note that screening 
would only reduce the contribution from this effect). 
We find 
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where the first equation suggests that the mean free 
path I for surface roughness scattering is typically 
larger than the sample size. These results, natu- 
rally, strongly depend on how one models the surface 
roughness. Recently, an alternate model proposed by 
Ref. ((17I ) suggested that in some special circumstances 
ripples can mimic Coulomb scatterers, although their 
phonon model additionally predicts strong tempera- 
ture dependence (above a certain quenching temper- 
ature of about 100K) - an effect that has not been 
observed in the experiments. We do not believe that 
the ripples in graphene are playing an important direct 
quantitative role in determining the carrier mobility in 
currently available samples, but rather scattering due 
to unintentional charged impurities is the dominant 
scattering mechanism in graphene. 

In Fig. 1 we show the effect of changing substrate 
dielectric constant on graphene transport proper- 
ties. The blue circles show the experimental data of 
Ref. (2) for a high mobility sample, and the red solid 
curve (k = 3.9) uses the experimentally reported val- 
ues in that paper for the undetermined parameters 
in our theory. A single fit parameter is still neces- 
sary to obtain the overall scale of the short-range 
impurities. Shown also are theoretical predictions for 
different dielectric constants. For k < 5, Coulomb 
scattering dominates (as is the case in most current 
experimental samples), while for larger dielectric con- 
stant, Coulomb scattering is more strongly screened 
making short-range scattering more dominant (see 
broken curve with k = 80). For intermediate values 
of dielectric constant, we have a crossover from low- 
density Coulomb-scattering-dominated behavior to 
high-density short-range-scattering dominated (see 
dashed curve with re = 10). 



3. Discussion and Conclusion 

We have examined three different sources of scat- 
tering in graphene and treated them in the RPA- 
Boltzmann regime which we believe provides the basic 
formalism to understand current graphene experi- 
ments. This formalism strongly indicates that close to 
the Dirac point, the transport properties are domi- 
nated by Coulomb impurity scattering. 

The basic picture is quite intuitive. As described 
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Fig. 1. Figure showing the dependence of conductivity on 
substrate dielectric constant ft. Filled blue circles show ex- 
perimental data of Ref. (Q). Solid, dashed and broken lines 
show theoretical results for ft = 3.9, 10, 80 respectively. For 
large ft short-range scattering dominates while for small ft 
Coulomb scattering dominates. The inset shows that for a 
fixed gate voltage, the conductivity has a non-monotonic 
dependence on ft which is a consequence of the competition 
between short-range and long-range scattering. 



above, within the RPA-Boltzmann formalism, charged 
impurities in the substrate provide a contribution to 
the conductivity that is proportional to the number of 
carriers in the graphene sheet. Additionally, these same 
charged impurities generate potential fluctuations that 
are screened by any carriers in the graphene sheet. A 
self-consistent procedure @) is used to determine the 
residual carriers that contribute to Boltzmann conduc- 
tivity even as the external gate voltage is tuned through 
the duality point where the majority of carriers changes 
from electrons to holes. This mean-field residual car- 
rier density n* provides for a minimum conductivity 
plateau and is in quantitative agreement with recent 
experimental data (QEI). In this context, the current 
work demonstrates that other scattering mechanisms 
do not alter this picture. Additionally, this work makes 
several predictions that can be tested experimentally if 
one is able to change the substrate onto which graphene 
has been cleaved. Experiments have already been per- 
formed on graphene sandwiched between various kinds 
of plastics using a transfer printing method (|l9T ) and 
experiments in the near future will certainly test the 
predictions we make in this work. 
This work is supported by US-ONR. 
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